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1.  INTRODUCTION 

Detecting  the  number  of  signals  and  estimating  the  parameters  of  the  damped  exponential 
signals  are  important  problems  in  signal  processing.  We  formulate  the  problem  as  follows. 
Let  yi ,  V2,  ■  •  •  >  Vn  be  a  sample  of  size  n,  where  yt  is  given  by 

yt  —  ^  ockexp(—skt  +  ttnfkt)  +  et  (!•!) 

k=i 

Here  at's  are  unknown  complex  numbers  called  the  amplitude  of  the  k‘h  signal,  A’s  are 
distinct  real  numbers  lying  between  0  and  1,  .*■; s  are  the  damping  factors  and  are  positive 
real  numbers,  i  =  t/=T.  {a.}  la  a  sequence  of  independent  identically  distnbuted  random 
variables  with  mean  zero  and  finite  variance  for  both  the  real  and  the  imaginary  part.  'Hie 
real  and  imaginary  part  of  {r,}  are  assumed  to  be  independent  and  normally  distnbut^ 
M,  the  number  of  signals  is  also  assumed  to  be  unknown.  Given  the  sample  of  size  n,  the 
problem  is  to  estimate  the  unknown  parameters  ak,  sk,  fk  for  k  —  1,  . .  •  M  and  M  also. 

The  estimation  of  the  parameters  of  a  damped  exponential  model  (1.1)  is  an  old  problem 
(Kay  1987)  and  the  readers  are  referred  to  Stoica  (1993)  for  an  extensive  list  of  references. 
Lots  of  methods  for  estimating  the  frequencies  have  been  proposed  by  researchersover  e 
last  twenty  years.  Among  the  notables,  are  the  methods  of  Ernkson  etaL{199^Kv 
(1984)  Kundu  and  Mitra  (1995),  Stoica  and  Nehorai  (1989),  Stoica  et  al.  (1989),  Tuft 
Kumaresan  (1982)  and  Yan  and  Bressler  (1993).  All  these  methods  of  estimation  assume 
that  the  number  of  signals  M  is  known.  The  aim  of  this  paper  is  to  estimate  the  nmnber  of 
signals  M,  which  is  usually  unknown,  under  the  assumption  that  the  number  of  signals  can 

be  at  most  K,  which  is  known  in  advance. 

Wax  and  Kailath  (1985)  developed  information  theoretic  criteria  for  detecting  the  num¬ 
ber  of  signals  received  by  a  sensor  array.  Fuchs  (1988)  developed  a  criterion,  based  on  the 
perturbation  analysis  of  the  data  auto  correlation  matrix,  for  detecting  the  nunJbe^  °  ^ 
soids.  More  recently  Reddy  and  Biradar  (1993),  following  the  information  theoretic  approach 
to  model  selection  developed  a  criterion  for  detecting  the  number  of. da™pe?/u*da™^  * 
ponentials.  The  detection  performance  of  these  criteria  were  compared 
(1988)  and  their  results  showed  that  the  Minimum  Description  Length  (MDL)  criterion  as 
developed  by  them  performs  nearly  same  as  that  of  Fuchs  (1988).  A  more  general  infor¬ 
mation  theoretic  criterion  in  model  selection  has  been  proposed  by  Zhao, Rested 
Bai  (1986a,  1986b)  called  the  Efficient  Detection  Criterion  (EDC).  Rao  (1988)  suggested 
to  use  EDC  to  estimate  the  number  of  signals  for  damped  or  undamped  case  but  he  did 
not  perform  any  numerical  experiments.  It  is  known  (Bai  et.  al;  1987)  that  the  EDf  g' 
consistent  estimates  for  estimating  the  number  of  signals  in  undamped  exPonantia 
although  the  same  result  is  not  applicable  for  damped  exponential  model.  Kundu  (1992) 
gave  a  detailed  comparison  of  the  different  information  theoretic  criteria  for  estimating  he 
number  of  undamped  signals,  but  nowhere  at  least  not  known  to  the  authors,  the  comparison 
of  the  different  information  theoretic  criteria  exist  for  damped  exponential  model. 
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Note  that  for  the  damped  exponential  model  the  data  sequence  is  pure  noise  as  the  sample 
size  goes  to  infinity.  Therefore,  one  can’t  obtain  any  asymptotically  consistent  estimate  of 
the  number  of  signals.  However,  when  the  damping  factor  is  not  that  first,  it  is  hoped  that 
some  good  detection  criterion  can  surely  be  obtained  by  suitable  algorithms,  which  should 
be  able  to  estimate  the  number  of  signals  reasonably  well.  That  is  the  main  aim  of  this 

paper. 

For  the  undamped  exponential  models  all  the  information  theoretic  criterion  can  be  writ¬ 
ten  in  the  form  (2.8),  where  Cn  represents  penalty  function.  It  has  to  satisfy  the  conditions 
given  in  (2.7).  Note  that,  the  penalty  function  Cn  goes  to  infinity  for  the  undamped  model 
to  give  consistent  estimate  of  the  number  of  signals.  For  the  damped  model  if  Cn  goes  to 
infinity,  then  for  large  sample  size  any  criterion  will  underestimate  the  number  of  signals. 
In  fact,  the  penalty  function  should  go  to  zero  as  n  tends  to  infinity.  We  modify  Cn  for  the 
damped  model  and  propose  the  modified  information  theoretic  criteria  where  the  penalty 
function  depends  on  the  amplitude  as  well  as  the  damping  factor.  If  there  is  no  damping  fac¬ 
tor  it  coincides  with  the  information  theoretic  criteria  for  the  undamped  model.  We  obtain 
the  probability  of  the  wrong  detection.  The  probability  of  wrong  detection  depends  on  the 
unknown  parameters.  We  propose  to  use  the  bootstrap  techniques  to  estimate  the  probabil¬ 
ity  of  wrong  detection  for  a  particular  penalty  function.  Once  we  estimate  the  probability  of 
wrong  detection,  we  choose  that  penalty  function  for  which  the  wrong  detection  is  minimum. 
Some  simulations  are  performed  to  see  the  effectiveness  of  the  proposed  criterion. 

The  organization  of  the  rest  of  the  paper  is  as  follows.  In  section  2  we  introduce  different 
information  theoretic  criteria  and  propose  the  modified  efficient  detection  criteria  for  the 
damped  model.  The  practical  implementation  procedures  are  provided  in  Section  3.  In 
Section  4  we  present  the  numerical  experiments  and  finally  we  draw  conclusions  in  Section 

5. 

2.  DIFFERENT  INFORMATION  THEORETIC  CRITERIA 

In  this  section  we  discuss  the  different  information  theoretic  criteria  for  estimating  the 
number  of  signals  of  the  damped  exponential  signal  models.  We  introduce  the  Akaike  Infor¬ 
mation  Criteria  (AIC),  Minimum  Description  Length  (MDL)  criteria  and  Efficient  Detection 

Criteria  (EDC). 

Let  t/i, . . . ,  yn  be  a  sample  of  size  n  from  the  model  (1.1).  Let  P,  be  the  parameter  that 
ranges  over  all  possible  number  of  signals,  i.e.  P  6  {1,...,A}.  Then  the  joint  density 
function  of  the  data  set  can  be  written  as 

f(y\°p)  =  ^^exp  (“|  £ |yt "  ^ 9p )|2)  (2,1) 

where 

dp  =  (ai,...,  ap,si,...,sp,fi,...1fP) 
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and 


»  r 

Ht{Op)  =  J2  Oikexp(-skt  +  i2ixfkt) 

k=i 

We  now  formulate  the  problem  as  follows;  Given  a  set  of  n  observations  and  a  family  of 
models  {f{y\0p)]P  =  1, that  is  a  parameterized  family  of  probability  densities 

f(y\9p),  our  problem  is  to  select  the  true  one. 

Posed  this  way,  this  problem  is  perfectly  suited  for  using  different  information  theoretic 
criteria  such  as  AIC,  MDL  or  the  EDC.  The  AIC,  MDL  and  EDC  criteria  are  known  as 
penalized  likelihood  method  in  the  general  statistical  literature.  Here  a  penalty  function 
is  subtracted  from  the  log-likelihood  before  it  is  maximized.  This  serves  to  penalize  or 
discourage  the  addition  of  more  and  more  parameters.  In  this  set  up  the  best  model  would 
be  one  for  which  the  penalized  likelihood  is  maximum.  For  the  general  problem  on  this  topic 
one  can  refer  to  Akaike  (1973,  1974,  1978),  Hannan  and  Qumn  (1979),  Rissanen  (1978), 
Schwartz  (1983)  arid  Zhao  et  al.  (1986a,  1986b). 

Akaike  (1973, 1974)  proposed  the  Akaike  Information  Criterion  (AIC).  The  AIC  suggests 
choosing  M,  an  estimator  of  M,  which  minimizes  the  following  expression; 

AIC{P)  =  - logf{y\9p )  +  d{9P)\  (2*2) 

for  P  =  where  9P  is  the  maximum  likelihood  estimator  (MLE)  of  9P  and  d(9P)  is 

the  number  of  independent  parameters  of  the  parameter  vector  9P. 

Akaike’s  basic  idea  was  to  choose  the  model  that  minimizes  the  mean  of  the  Kullback- 
Leibler  distance  between  the  true  density  f(y\9P)  and  the  estimated  density  f(y\9P).  Smce 
the  distance  is  unknown,  he  proposed  to  estimate  it  by  the  log-likelihood  of  the  MLE.  The 
second  term  in  (2.2)  was  added  to  make  the  log-likelihood  at  the  MLE  an  unbiased  estimator 
of  the  Kullback-Leibler  distance. 

In  the  exponential  signals  model,  with  the  assumption  of  the  Gaussian  error  the  AIC 

takes  the  following  form;  .  . 

AIC(P)  =  —nlogRp  -  8 P  (2-3) 

(see  Rao;  1988),  where  RP,  denotes  the  minimum  value  of 

EIk-MMI2  (2-4) 

t=  1 

and  the  minimization  is  performed  with  respect  to  ai, . .  • ,  aP,  s\, . . . ,  sP,  /i,  •  •  • ,  fp- 


Minimum  Description  Length  (MDL)  criterion  was  introduced  by  Rissanen  (1978)  The 
basic  idea  is  that  the  best  model  is  the  one  that  provides  the  shortest  description  of  the 
data.  It  has  been  shown  (Rissanen;  1983)  that  for  large  samples  this  criterion  leads  to  the 

selection  of  the  model  that  minimizes 

MDL(P)  =  - logf{y\9P )  +  i d{9P)logn  (2-5) 


4 


fox  P  =  l,...,K,  where  f(y\9P)  and  d(9P)  are  as  defined  before. 

Schwartz  (1978)  suggested  a  model  selection  approach  based  on  Bayesian  arguments. 
Assuming  a  priori  probabilities  for  every  competing  model,  he  proposed  selecting  the  model 
that  maximizes  the  posterior  probability.  It  has  been  shown  that  for  a  model  belonging 
exponential  family,  the  maximization  of  the  posterior  probability  leads  to  the  minimization 

of  the  criterion  given  by  (2.5)  asymptotically. 

The  Efficient  Detection  Criterion  (EDC)  method  of  Zhao  et  al.  (1986a,  1986b)  consists 
of  choosing  as  an  estimator  of  M,  the  number  M,  which  minimizes 

EDC(P)  =  —logf{y\0p)  +Cnd(6p)  (2*6) 


for  P  =  1, . . . ,  K,  where  C„’ s  are  such  that 

Q 

lining  oo—  =  0  and  limT 


3  loglogn 


In  the  exponential  signal  model,  with  the  assumption  of  Gaussian  error,  the  EDC  takes  the 
following  form  (see  Rao;  1988) 

EDC(P)  =  -nlogRp  -  C„(8P)  (2-8) 

Observe  that  MDL  criterion  is  a  special  case  of  the  EDC  For  MDL,  C„  takes  the  value  | 
loe  n  in  (2  8).  The  estimators  of  M,  obtained  from  (2.8)  are  strongly  consistent 
undamped  exponential  model.  For  a  detailed  proof  of  the  consistency  for  the  undamped 
model  see  Bai  et  al.  (1987).  For  the  damped  model  however  the  consistency  results  do  not 
hold,  therefore  it  is  important  to  observe  the  behavior  of  the  different  information  theoretic 
criteria  in  this  situation  at  least  for  small  samples. 

Now  we  try  to  analyze  what  kind  of  problem  we  might  encounter  if  we 
for  estimating  the  number  of  signals  for  damped  exponential  model.  Note  that,  (2.7)  P 
Cn  tends  to  infinity  as  n  tends  to  infinity.  Suppose,  M  is  the  correct  order  model,  then  Gn 

should  be  such  that 


EDC{M)  <  EDC(P)]  for  P  =  1,  •  •  • ,  K,  P  ±  M. 
Now  (2.9)  implies 

nlogRM  +  Cn(8M)  <  nlogRP  +  Cn(8P);  for  P  =  1,  •  •  • , K, P 
Since  Ri  >  R2  >  ■  •  •  >  &k  almost  surely,  (2.10)  implies  that  C„  must  satisfy 


p  =  ±  M.  (2.10) 


<  8Cn  <  nlog 


Rm-i 

Rm 


(2.11) 
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For  undamped  model 


,.  Rm 
lim  — - 

n-*°°  Rm+1 

=  1, 

and 

Rm-  1  ^  , 
lim  — —  >  1 

n->°°  Rm 

(2.12) 

for  damped  model 

,.  Rm 
hm 

Rm+i 

=  1, 

and 

Rm-1 

lim  =1 

Rm 

(2.13) 

Because  of  the  damped  factor,  note  that  for  large  n, 


=  1  +  0  (e~5n)  (2.14) 

Rm  v  ' 

where  8  =  max  {si,...,sm}  >  0-  If  we  divide  by  loglogn  in  (2.11)  and  take  the  limit,  we 
obtain 

•  oo  =  8  lim  7 — 7— —  <  lim  nlog  (l  +  O  (e-<5n))  =  0  (2.15) 

n-too  loglogn  n~* 00  '  v  ' ' 

Therefore,  if  Cn  tends  to  infinity,  for  large  n,  (2.11)  may  not  satisfy.  On  the  other  hand  it 
looks  more  reasonable  that  the  penalty  function  should  be  more  if  the  amplitudes  are  more 
(suggested  by  a  referee).  Based  on  the  above  observations,  we  propose  the  following  modified 
EDC  (MEDC)  for  the  damped  model 

MEDC(P)  =  —nlogRp  -  ACne~6n{8P )  (2.16) 

here  A  =  max{Ai, . . . ,  AM)  and  5  =  max{8 u  . . . ,  8M}-  If  the  damping  factor  is  zero,  then 
MEDC  coincides  with  the  usual  EDC.  Note  that  we  need  to  know  A  and  8  to  implement 
MEDC  in  practice.  We  will  describe  that  in  the  next  section. 

3.  PRACTICAL  IMPLEMENTATION 
Consider  the  following  data  matrix, 

Vi  •  •  •  Vl 
A  =  : 

2/n-L+l  •  *  •  Vn 

Here  L  is  any  integer  such  that  K  <  L  <  N  —  K.  Let  us  denote  the  matrix  T  =  ^A*A, 
where  **’  denotes  the  conjugate  transpose  of  a  matrix  or  of  a  vector.  We  obtain  the  spectral 
decomposition  of  the  matrix  T  as  follows: 

L+ 1 

T  =  H  a,2u(u; 

t=l 

here  erf  >  ...  >  o\+l  are  the  ordered  eigenvalues  of  T  and  U,’s  are  the  normalized  eigenvalues 
corresponding  to  <5f . 
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Assuming  the  true  order  of  the  model  is  if  (the  maximum  one)^  we  “ 

K  damping  factors  and  the  K  amplitudes  say  <5i  >  . . •  >  &k  and  AU...,AK  p  ^ 
u  d<=in!  the  NSD  method  of  Kundu  and  Mitra  (1995)  from  T.  We  use  <5  =  <5i  and  A  — 
max  ",8  .  Zb  Note^hat  the  values  of  4  and  A  depend  on  L,  we  provide  some  suggestmns 

to  choose  L  in  the  next  section. 

For  a  given  choice  of  C„  and  from  the  estimated  A  and  4,  we  can  compute  MEDC(P) 
for  different  values  of  P  =  1 . K  and  choose  M  an  estunate  of  U  such  that  MEDC(M) 

is  minimum. 

Note  that  we  have  a  wide  choice  of  C„,  but  we  would  like  to  choose  that  C„  so  that 
P(M  7^  M)  is  minimum.  First  let’s  compute  P(M  i=-  M). 


P(M  ±  M)  =  P(M  <  M)  +  P{M  >  M) 


M-l  JL  ->  . 

.  V  P{M  =  q)+  £  = 

~0  q=M+ 1 

P(MEDC(q)  -  MEDC(M)  <  0) 

9=° 

£  P(MEDC(q)  -  MEDC(M)  <  0) 

q=M+ 1 

P(nlogRq  —  nlogRM  >  ACne  8(M  —  q)) 

9=° 

'y'  P(nlogR\f  —  nlogRq  <  ACne  8 (q  —  M)) 

q=M+ 1 


(3.1) 


it  r  f„wp1v  P(M  ±  M)  depends  on  the  unknown  model  parameters.  Without  knowing  the 
orSr—  £  can’t  calculate  the  theoretical  probabilities.  We  would  like  to  estimate 
these  probabilities  with  the  help  of  the  given  sample  and  using  the  bootstrap  teduuqu  . 
idea  is  as  follows.  From  any  particular  realization  of  the  model ™  ^  by 

obtain  the ^cmrrKpondii^ei^emmlu^andeigenv^ctom^^Weesiumate^tlie^error^v^^  junct^ 

“estimate  the  order  of  the  model  as  Af(C  ).  We  generate  n  “^x  Ga— anjm 

variables  with  mean  zero  and  variance  a  ,  say  ci,.--,  n- 

sample  as  ,  ,  ,  „ 

yf  =  yt  +  6‘»  for  t  —  l,...,n. 

Assuming  M(Cn)  is  the  correct  order  model,  we  check  for  q  <  M(Cn),  whether 

nlog{Rq)  —  nlog(RM(cn))  >  ACne  8 (M(C„)  —  q ), 

or,  for  q  >  M(C„),  check  whether 

nlog(RM{cn))  ~  nlog(Rq)  <  ACne~5nS(q  -  M(Cn)\ 
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Repeating  the  process,  say  B  times,  we  can  estimate  (3.1).  Finally  we  choose  that  Cn  for 
which  the  estimated  P{M  ^  M)  is  minimum. 

Some  justifications  regarding  this  kind  of  bootstrap  estimates  of  (3.1)  can  be  given.  Note 
that  the  realization  of  yf  can  be  thought  of  coming  from  a  model  (1.1)  with  V(et)  «  2a2. 
Note  that  for  the  damped  exponential  model  -§*-  for  q  =  1, . . .  ,K  is  independent  of  a2. 

t  2  o  2 

Therefore  (3.1)  remains  invariant  if  we  change  the  error  variance  from  a  to  2a  . 

4.  NUMERICAL  EXPERIMENTS 

In  this  section  we  present  the  Monte  Carlo  simulations  done  for  small  samples  to  compare 
the  different  information  theoretic  criteria.  All  these  computations  have  been  done  on  HP  — 
9000,  machine  at  the  Indian  Institute  of  Technology,  Kanpur. 

We  consider  four  different  models  with  different  parameters  and  different  standard  devi¬ 
ations  of  the  error  random  variables.  The  four  models  are  given  as  follows; 

Model  1  yt  =  c*/4c(-0.0«-H2»(.5a)l)  +  e*/2e(-0.02t+i2n(A2)t )  +  €t 
Model  2  yt=  e'/4e(-0  0U+^(.52)t)  +  ^(-O.OZt+^.SO)*)  +  Ct 
Model  3  yt  =  e^V"00 "+&*{&)*)  +  e^e(-o.o3t+i2^A2)t)  + 

Model  4  yt  =  1+  e^e(-0.01t+<2.(.52)t)  +  c»/2e(-0.02M-i2,r(.50)t)  +  €t 

The  data  are  generated  using  the  different  standard  deviation,  viz  a  =  0.01,  0.1,  0.5  and  1.0 
and  with  different  sample  sizes  n  =  25,  50,  75  and  100.  The  random  deviates  are  generated 
with  the  help  of  the  IMSL  random  deviate  generator.  For  each  of  the  four  models  one 
hundred  replications  of  the  data  set  for  different  n  and  a  are  generated.  Observe  that  in 
Model  1  and  Model  2,  the  amplitudes  and  the  damping  factors  are  kept  fixed,  whereas  the 
difference  of  the  radian  frequencies  is  more  in  Model  1  than  in  Model  2.  Between  Model  1 
and  Model  3,  the  amplitudes  and  the  radian  frequencies  are  kept  fixed,  whereas  the  difference 
between  the  damping  factor  is  more  in  Model  3  than  in  Model  1.  Model  4  is  a  higher  order 
model  than  Models  1,2  or  3.  As  far  as  the  estimation  of  frequencies  are  concerned,  it  is 
known  (Kundu  and  Mitra;  1995)  that  it  is  difficult  to  estimate  the  parameters  in  Model  2 
than  in  Model  1  and  similarly  in  Model  1  than  in  Model  3.  No  such  comparison  can  be  made 
between  Model  2  and  Model  3.  Between  Model  4  and  Model  2  it  is  expected  that  Model  2 
will  be  easier  than  Model  4  as  the  number  of  parameters  are  more  in  Model  4  than  that  of 
Model  2.  It  is  expected  that  in  estimating  the  number  of  signals  also,  the  same  pattern  will 
exist. 

We  compare  the  usual  AIC  and  usual  MDL  with  the  proposed  MEDC.  Note  that  for  AIC 
and  MDL,  Cn  =  1  and  Cn  =  \logn  respectively  in  (2.8).  For  MEDC,  we  take  a  varied  choice 
of  Cn  satisfying  (2.7)  (except  when  Cn  =  1)  but  diverging  to  infinity  at  different  rates  from 
very  slow  to  very  fast.  The  different  choices  of  Cn  considered  are  as  follows:  Cn  =  1,  C„  =  tv 
,  Cn  =  n-5,  Cn  =  n-9  Cn  =  log  n,  Cn  =  ( logn )1,  Cn  =  ( logn)-5,Cn  =  ( logn )'9  Cn  =  (nlogn)1. 
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■c  -  (nloqn)  5  Cn  =  (nlogn)9  and  Cn  =  \  log  n.  It  is  assumed  that  for  all  the  four  models 
the  "maximum  number  of  signals  is  6.  Firs?  we  assume  that  the  -del  orier  is^  =  6^  Now 
using  the  modified  noise  space  decomposition  method  with  L  «  min 

estimates  of° A  and  5.  Using  that  a4  and  5,  we  compute  MEDC(P)  from  (2  16)  for  different 
values  of  F  =  1,  •••  6  for  a  particular  choice  of  Cn.  We  obtain  an  estimate  of  a  y  averaging 
the  last  (L  -  K)  eigenvalues  of  AHA  and  also  obtain  an  estimate  of  P{M  ^  M)  as  the 
method  suggested  in  the  previous  section.  We  take  B  =  100,  in  our  calculations  The  results 
are  reported  in  Tables  1-4.  We  report  the  percentage  of  under  estimate  (PUE^  P^entJ 
of  correct  estimate  (PCE)  and  the  percentage  of  over  estimate  (POE)  for  AIC,  MDL  and 
MEDC  over  five  hundred  replications. 

5.  CONCLUSIONS: 

In  this  paper  we  consider  the  problem  of  estimating  the  number  of  damped  exponential 
signals.  We  use  different  information  theoretic  criteria  for  estimating  the  number  of  sign  . 

We  consider  the  AIC,  MEDC  and  the  MDL  criteria  for  the  detection  problem.  It  is 
well  known  that  the  AIC  criteria  does  not  provide  the  consistent  estimates  in  general  mo  e 
selection  problem.  This  fact  is  weU  reflected  in  the  results  of  the  simul m  a 
1-4.  Comparing  the  Tables  1-4  it  is  observed  that,  although  the  MEDC  and  MDL  critena 
give  consistent  estimates  for  undamped  signals  model,  the  same  can  not  be  said  for  the 
damped  models.  It  is  well  known  (Wu;  1981  and  Kundu;  1994)  that  although  it P<«.Me 
to  estimate  consistently  the  parameters  of  the  undamped  exponentia1  model  but  it  ^  not 
possible  to  obtain  the  consistent  estimates  of  the  parameters  of  the  da  P  '  ... 

nm  be  surprising  if  we  look  at  the  damped  model  carefully.  From  the  model  it  is  dear  (if  the 
damping  fLtor  is  negative)  as  the  sample  size  n  increases  the  signal  component  ™mehes  to 
zero  and  we  left  with  the  error  components  only.  Therefore  even  if  we  increase  the  sa  p 
sizes  we  may  not  extract  any  more  information  about  the  signal  parameters  from  the  sample. 
In  fit  the  inconsistency  is  clearly  indicated  in  the  simulation  results.  It  is  observed  that 
the  number  of  correct  selections  by  different  methods  do  not  increase  for  a  fixed  a  as  n 
increases.  In  fact  for  AIC  and  MDL  in  many  cases  they  even  decrease.  For  fixed  n,  as  o 
decreases,  it  is  observed  that  for  MDL  and  MEDC,  the  performances  improve.  This  indicate 
the  consistency  of  the  MDL  and  MEDC  methods  as  <r  decreases  to  zero  for  fixed  n.  : It  is  ato 
observed  that  for  a  fixed  n  as  a  increases  the  methods  have  a  tendency  to  over  estimate  for 
models  1  and  3,  whereas  they  have  a  tendency  to  under  estimate  for  model  2  and  mo^4 
Comapring  the  tables,  it  is  observed  that  in  most  of  the  cases,  for  fixed  n  , 
of  correct  detection  is  more  in  Table  1  than  in  Table  2,  which  is  not  very  surprnnng j»  *• 
difference  of  the  radian  frequencies  (|/i  -  h\)  is  more  in  model  1  than  in  model  2.  This 
was  also  seen  for  the  undamped  signals  by  Kundu  (1992).  Interestingly  the  °f 

detection  in  Table  3  is  more  or  less  same  as  that  of  Table  1  although  the  difference  in  the 
damping  factor  (Is,  -  s2|)  is  more  in  model  3  than  in  model  1.  The  difference  in  performance 
is  more  marked  if  the  two  models  differ  significantly  with  respect  to  the  frequencies.  The 
performance  of  most  of  the  methods  is  much  better  for  model  1  than  that  of  model  2  if  a 
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.  n  m  Between  Model  2  and  Model  4,  the  behavior  are  quite  similar  in  nature  for  all 
^tXhe“o»ied,  although  the  number  of  correct  detection  .  more  »  Mode!  2 

compared  to  Model  4. 

Now  comparing  the  three  methods  it 

particular  model.  Our  simulations  s  ow  prominent.  MDL  criterion  works 

never  exceeds  .45  also  the  inconsistency  of  the  y  P  ...nance  is  ygt,  and  the 

reasonably  well  if  the  error  variance  is  not  very  h1^  H  the^  performance 

difference  of  the  radian  frequencies  is  small  (Model *  “d  »dode 1  4)  thjn  P 
of  MDL  is  also  very  poor.  MEDC  work  very  well  if  the  rror  var  an  is  low^ 

i  /-kef  on  tn  inn  nercent  for  all  the  models  considered  if  o  —  T 

almost  90  to  1UU  percent  lor  a  frequencies  are  close  to  each  other.  It 

high  the  performance  drops  significantly  if  the  radian  irequenue. 

Although,  eventually  as  n  tends  to  infinity  MEDC  also  will  giv  ^  knQwn  methods. 
least  for  finite  sample  it  works  reasonably  well  and  be  a  jp  MDL 
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